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ABSTRACT 

We have analysed the XMM-Newton spectra of SS433 using a standard model of 
adiabatically and radiatively cooling X-ray jets. The multi-temperature thermal jet 
model reproduces well the strongest observed emission line fluxes. Fitting the He- and 
H-like iron line fluxes, we find that the visible blue jet base temperature is » 17 keV, 
the jet kinetic luminosity Lj, ~ 2 • 10 39 erg/s and the absorbing column density Njj ~ 
1.5 • 10 22 cm~ 2 . All these parameters are in line with the previous studies. The thermal 
model alone can not reproduce the continuum radiation in the XMM spectral range, 
the fluorescent iron line and some broad spectral features. Using the thermal jet-plus- 
reflection model, we find a notable contribution of ionized reflection to the spectrum in 
the energy range from ~ 3 to 12 keV. The reflecting surface is highly ionized (£ ~ 300), 
the illuminating radiation photon index changes from the flat spectrum (T m 2) in the 
7-12 keV range to T w 1.6 in the range of 4-7 keV, and to T < 1 in the range of 
2-4 keV. We conclude that the reflected spectrum is an evidence of the supercritical 
disc funnel, where the illuminating radiation comes from deeper funnel regions, to be 
further reflected in the outer visible funnel walls (r > 2 ■ 10 11 cm). In the multiple 
scatterings in the funnel, the harder radiation > 7 keV may survive absorption, but 
softer radiation is absorbed, making the illuminating spectrum curved. We have not 
found any evidences of reflection in the soft 0.8-2 keV energy range, instead, a soft 
excess is detected, that does not depend on the thermal jet model details. However the 
soft component spectrum is basically unknown. This soft component might prove to be 
the direct radiation of the visible funnel wall. It is represented here either as black body 
radiation with the temperature of 9u> « 0.1 keV and luminosity of Lw ~ 3 • 10 37 erg/s, 
or with a multicolour funnel (MCF) model. The soft spectral component has about 
the same parameters as those found in ULXs. 

Key words: massive close binaries - stellar mass black holes: individual: SS 433 - 
X-ray spectra: supercritical accretion discs - ULXs. 



1 INTRODUCTION 

SS 433 is the only known persistent superaccretor in the 
Galaxy - a source of relativistic jets l|Fabrikall2004l . for re- 
view) . This is a massive close bina ry, where the compact star 
is m os t probably a black hole (IGies. Huang, fe McSwain 



2002 ; ICherepashchuk et all 120051; iHillwig fe Giesl 1200 



20081 ). 



SS433 

-,40 , 



Blundell. Bowler, fe Schmidtobreickl 

trinsic luminosity is estimated to be ~ 10 40 erg/s, 
with its maximum located in non-observed UV re- 
gion. In effect this is a very blue, but heavily absorbed 
object; its estimated te mperature depends on the ac- 



cretion disc orientation dMurdin. Clark. fe Martini fl98tj|: 



Cherep ashchuk. ' Aslanov. fe Kornilovl 19821 ; Dolan et al.1 
19971 ). being ~ 50000 - 70000 K when the disc is the most 



open to the observer. It is important that almost all the 
observed radiation is formed in the supercritical accretion 
disc, and the donor star contributes less than 20 % of the 
optical radiation. The system's extreme luminosity suggests 
that the black hole's mass is ~ 10 Mq . Such a big energy 
budget of the object is supported by a very well measured 
kinetic luminosity of the jets, ~ 10 39 erg/s, both in direct 
X-ray and optical studies of t he jets and in the studies of 
the jet-powered nebula W 50 (|Fabrikall200i ). At the same 
time we know that practically all the energy at the accretion 
onto a relativistic star is released in X-rays. This means 
that the observed radiation of SS 433 was thermalised in 
the strong wind coming from the supercritical disc. 

The wind from a supercriti ca l dis c was first 
predicted by IShakura fe Svunvaevl (|l973l ) and later 
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confirmed in radiati on -hydr ody namic sim ul ations 
llEggum, Coroniti. fc Kat3 Il988l ; lOkuda et af] 120051 ; 
lOhsuga et alT feoOS). These ideas led to a prediction that 
SS 433, being observed face-on appears as an extremely 
bright X-ray source, and we may expect an appearance 
of a new type of X-ray sour c es in galaxies (|Kat d 1 19871 : 
iFabrika fc Meschervakovl |2000| . 1200 ll ) - the face-on SS 433 
stars. It is quite possible that the new type of X - ray so urces, 
ULXs (ultraluminous X-ray sources, iRobertsI (|2007h ^ arc 
SS 433-like objects observed nearly face-on. Their observed 
X-ray luminosities are 10 39 - 3 • 10 41 erg/s and they arc 
certainly related to the massive star population. The disc 
orientation in SS433 is about edge-on, the angle between 
the disc plane and the line of sight is 10 ± 20° , due to the 
precessional variations. Therefore we have no chance to 
observe the funnel in the supercritical disc directly. 

The supercritical disc radiation is not isotropic. If one 
takes into account the geometrical beaming in the disc 
funnel with the (h alf) opening angle i?/ ~ 25° — 30° 
jOhsuga et al.| [2005'l and t hat a supercritical disc surfa ce ra- 
diatio n is local Eddington l|Shakura fc Svunvaevlll973T ). one 
finds (I Fabrika et al .112006! ) . that the observed X-ray luminos- 
ity of ~ 10 41 erg/s is e xpected for the face-on SS433 star . 
The recent data sh ow (|Stobbart. Roberts, fc Wilms! l200rj : 
Bergh ea et al.ll2008l ) that ULXs posses curvature and rather 
flat X-ray spectra, which are difficult to interpret with a 
single-component or any other simple model. The supercrit- 
ical accretion di scs are expected to h ave flat (yF u cx v°) 
X-ray spectrum l|Poutanen et al.ll2007h . because the energy 
release Q(r) in the discs must be Q(r) cx r~ 2 (T(r) cx r -1 / 2 ) 
to make the disc thick due to the radiation pressure. 

The X-ray lu minosity of SS433 is ~ 10 36 erg/s (at a 
distance of 5.5 kpc l|Blundell fc Bow ler 2004)), four orders of 
magnitude less than the bolometric luminosity. It is believed 
that all the X-rays come from the cooling X-ray jets. The jet 
may be easily accel erated by the radiat ion in the hydrody- 
namic funnel (e.g., lOhsuga et al.l l2005) of the supercritical 
disc. The jet velocity value, Vj ~ 0.26c, and its unique sta- 
bility, where the velocity does not depend on the activity 
state, indicate that the jet accel eration must be controlled 
by the line-locking mechanism l|Shapiro. Milgrom. fc Reesl 
ll986l : lFabrikall2004l ). The funnel has to be relatively trans- 
parent for the accelerating radiation, as it was confirmed 
in rad iation-hydrodynamic simulatuons (e.g., lOhsuga et al.l 
120051 ). 

These reasonings, however, bear a problem, why do 
not we observe any funnel radiation, missing four orders 
of magnitude in the X-rays? Even a small amount of gas 
in the most outer part of the funnel, where we observe 
the X-ray jet base, may reflect and scatter some part 
of the direct ~ 10 40 erg/s of the funnel radiation. The 
structure of the funnel and the jet acceleration/collimation 
region is unknown. Apparently, the direct funnel radi- 
ation is entirel y blocked for the observer. The 'stan - 
dard' jet model dBrinkmann et aL 19881: Kotani et al.ll 199 



IMarshall. Canizares. fc Schul j 120021 : iFilippova et aljBoO 
suggests that all the observed X-ray radiation of SS 433 is 
formed in the cooling X-ray j ets. However, an analysis of the 
latest XMM observations (jBrinkmann. Kotani. fc Kawail 
120051 ) led the authors to a conclusion that the 'standard' 
jet model can not produce the observed X-ray continuum. 
In the 'standard' jet model there are two conical anti- 



parallel jets, considered identical. The jet's gas is optically 
thin, being in collisional ionization equilibrium and cooling 
adiabatically (or adiabatically plus radiatively) . The jets are 
observed beginning from the jet base tq (the distance be- 
tween the base of the visible jet and the black hole, different 
for the red and blue jets), which depends on the preces- 
sional phase and on the eclipses by the donor star. The gas 
temperature at the jet base To(ro) is estimated from ob- 
servations. Precessional a nd orbital variability of SS433 is 
well-known (|Fab rika 2004|), the jets (and both the accretion 
disc, and the accretion disc wind) precess with a 164-day 
period and an amplitude of ±20° . At the precessional phase 
ip ~ the disc is the most open to the observer, the angle 
between the jets and line of sight is ~ 60° (we refer here 
to the approximate values because there are nutation-like 
and sporadic variabilities, both of about 3 — 5°). The ac- 
cretion d isc rim (or the opaqu e part of the wind) is thick 
h/r ~ 1 l|Filippova et al.ll2006h . Comparing the amplitudes 
of the pre cessional and orbital variab ilities in different X- 
ray bands, ICherepashchuk et al.l (|2005l ) found that both the 
donor star and the outer disc rim have about the same size, 
Therefore, one may expect the same radius for the approach- 
ing jet base ro ~ 10 12 cm. 

In X-ray obse r vations with ASCA and Chandr a 
dKotani et al. I Il996l : IMarshall. Canizares. fc Schulzl 120021 ). 
dozens of X-ray emission lines formed both in blue and 
red jets were resolved. The lines are variable both in po- 
sitions (in accordance with the jet kinematic model) and 
in intensities, the strongest are He- and H-like Fe XXV Ka 
and Fe XXVI La iron lines. This allowed the 'iron line di- 
agnostics' of the gas temperature at the jet bases by mea- 
suring the line ratio Fe XXV/FeXXVI. The temperature at 
the base of the jet was es t imated as 8p = kTp = 10 -j- 
20 keV dKotani et al] Il996l : IMarshall. Canizares. fc Schulzl 
l2002l : lNamiki et al.ll2003l). In the latest XMM- Newton obser- 
vations lBrinkmann. Kotani. fc Kawail (|2005l ) estimate the 
temperature of about Go ~ 17 ±2 keV, with the main uncer- 
tainty coming from the form of the underlying continuum. 
Fitting the cont inuum with a thermal b remsstrahlung model 
(GINGA data, iBrinkmann et al.lll99ll ) gives a notably big- 
ger value, 9o > 30 keV. The jet base temperature determined 
both in the line diagnostics and in the continuum fits, drops 
notably in eclipses, what confirms th e cooling jet model . 
Analysing the numerous ASCA data, iKotani et al] (|l996l ) 
found that the farther part of the receding jet is weakened 
(probably absorbed in the external gas located in the disc 
plane) and Nickel is highly overabundant (~ 10 times) in the 
jets. The l ast finding has been confirmed in th e XMM ob- 
servations l|Brinkmann. Kotani. fc Kawail [20051 ) . Note that 
some discrepancies may exist as the system SS 433 is highly 
variable and its appearance depends strongly on the preces- 
sional phase. 

IBrinkmann. Kotani. fc Kawail (2005) found that the 
XMM SS 433 spectrum can not be fitted with any sim- 
ple continuum law. The numerous lines point at the ther- 
mal model, but a strong continuum curvature indicates the 
Comptonization. The best formal model of the additional 
continuum component is a broken power law with a break 
at ~ 7.1 keV, where the low energy PL rises with energy 
(r ~ —1), the high energy part is rather steep (T ~ 4), how- 
ever, the parameters are not usually well determined in the 
fits. The overall continuum is too hard to be bremsstrahlung, 
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the electron-electron bremsstrahlung can not fit it even 
for the temperatures of 40 keV. A highly absorbed addi- 
tional thermal component giv e s a to o strong ~ 7 kev edge. 
iBrinkmann. Kotani. fc Kawail |2005h discuss a model with 
a single, very broad ~ 7 kev line, which may be formed 
in the innermost hot jet (otherwise hidden) and Compton 
down-scattered by colder material. These uncertainties in 
continuum make problems for the iron-line diagnostics of 
the jet base temperature. 

Practically in all the studies of SS 433 X-ray spectra the 
jet mass loss Mj and the jet kinetic luminosity Lk = MjVj /2 
were determined. To estimate Lk one needs to adopt the jet 
(half) opening angle $j. The mass loss rate derived from 
the X-ray line intensities depends strongly on the open- 
ing angle adopted, i. e. gas density at the jet base, be- 
cause the emission measure is oc n 2 . The kinetic luminos- 
ity values derived from X-ray spectra were quite diverse, up 
to 10 42 erg/s. However, in two latest studies, one from the 
Chandra data (|Marshall. Canizares. fc Schulz 2002) a nd the 
other from the XMM data (Brinkma nn. Kotani. fc Kawail 



l2005ft , the values are ~ 3 • 10 38 and ~ 5 ■ 10 39 erg/s re- 
spectively. They are even closer to one another, match- 
ing different distances to S S433 adopted by these authors. 
iPanferov fc Fabrikal (|l997ft found Lk ~ 10 39 erg/s from the 
jet Balmer line intensities in op tical spectra. The mass 
loss rate in the jets was estima t ed (IZealev. Dopita. fc Malinl 
Il98d : iFabrika fc Borisovl 1 19871 ; iDubner et all Il998h on "the 
base of W 50 nebula powered by the jets with approximately 
the same result, Lk ~ 10 39 erg/s. Thus the jet kinetic lumi- 
nosity was measured with relatively good accuracy. 

The jet opening angle was directly measured in th e 
studies of X-ray and optical l ine w idths (|Fabrikal |2004| ) . 
iMarshall. Canizares. fc Schulzl (|2002T | found &j = 0.61 ± 
0.03° from the Chandra spectra. They assumed that the den- 
sity is uniform through the jet c one's cross section a nd i?, 
is the jet's cross sectional radius. iBorisov fc Fabrikal (|l987l i 
found -dj — 0.82 ±0.14° from optical spectra, assuming that 
the emitting gas density is Gaussian-distributed in the cone's 
cross secti on. We give here the values o f i?, r ecalculated, 
using the [Marshall. Canizares. fc Schulzl (2002) definition. 
Such a coincidence in the X-ray and optical jet opening an- 
gles makes the jets pure ballistic from X-ray (~ 10 12 cm) 
to optical (~ 10 cm) regions. In optical spectra the jet 
opening angle has been estimated from Ha line profiles dur- 
ing 8-day long observations, the nutation broading has been 
taken into account and possible jet sporadic activity was 
diminished. In other C handra observations (jNamiki et al.l 
120031 ; lLopez et alj|2006ft the opening angle has been found 
twice as big and it was greater in the higher temperature 
(Fe), than in the lower temperature (Si) parts of the jets. 
Regarding comparatively short time scales of the X-ray ob- 
servati ons and probable li ne broading due to electron scat- 
tering (jNamiki et alj|2003l ) in the hottest parts of the jet, we 
may adopt -dj m 0.7° as a good estimate of the jet opening 
angle. 

In this paper we analyse X-ray spectra of SS 433 using 
XMM data and the 'standard' jet model. We do not try 
to determine the jet kinetic luminosity, because the X-ray 
continuum is com plex and can not be explained in the 'stan- 
dard' jet model l|Brinkmann. Kotani. fc Kawail [20051 ) . In- 
stead we adopt these two main parameters - Lk = 10 39 erg/s 
(Mj = 3.3 • 10 19 g/s) and 2<&j = 1.5°, as well established and 



known. We find that the line fluxes can be well matched 
with the 'standard' jet model and study the additional com- 
ponents of the X-ray spectrum of SS 433. 



2 OBSERVATIONS AND DATA REDUCTION 

We used a public archive of the XMM-Newton Observatory^] 
There are 8 observations of SS433 in the archive, where 
the target was in the detector's centre. They were carried 
out in April 2003 (4 observations at the precessional phase 
ip ~ 0.8), October 2003 (2 observations at ip w 0), and two 
single observations in March 2004 and April 2004. In this 
paper we concentrate on the study of X-ray jets and the thick 
accretion disc (its wind). Then we have selected only those 
observations, where disc was not eclipsed by the donor and 
was the most open to the observer. This is the observation 
taken on 19 October 2003 in orbit 707 (ip = 0.04) in the 
orb ital phas e <£ = 0.28 (for orbital and precessional phases 
see iFabrikal (|2004h 1. For comparison purposes we selected 
the observation from 6 April 2003 in orbit 609, taken at the 
intermediate inclination (ip ~ 0.84) of the disc and at about 
the same orbital phase (p = 0.34. The last observation is 
relatively short, hence we added the following observation 
taken in the same accretion disc precessional position, in 
orbit 610 (ip = 0.49) for a control check. 

These data allow us to study the jets and the disc wind 
at two different orientations of the disc/jets and at about 
the same orbital phase. With such a strong mass transfer in 
the system, the jet and the accretion disc funnel visibility 
behind the wind photosphere (the outer disc rim) may be 
dependent on the orbital phase. 

We use the EPIC PN observations to get the highest 
signal-to-noise data. During all the observations the EPIC 
PN camera was operated in the Small Window mode with 
a medium filter. The PN data were reprocessed using the 
XMMSAS version 6.1. The actual Life Time of the detector 
after a correction for Good Time Intervals was 11.4 ksec in 
orbit 707, 4.4 ksec in orbit 609 and 5.6 ksec in orbit 610. 



3 JET MODEL 

To model the jet spectrum we use the standard approach 
which was used by many authors t o analyse the X-ray spec- 
tra of SS 433 dKotani et all Il996t IBrinkmann et aL Il99ll ; 



IBrinkmann. Kotani. fc Kawail 120051 ; iFilippova et alj 12006ft . 
This standard multi-temperature model describes well both 
the jet line spectrum and its variability. The jet gas moves 
in ballistic trajectories with a velocity v = 0.26c. In such a 
case the gas concentration can be written as follows: 



n(r) 



no 



where no 



AL 



(1) 



Here ro is the visible jet base (the distance between the base 
and the black hole), fl,- = 2ti"(1 — cos^j) is the solid opening 
angle of the jet, Mj - the mass loss rate in the jet, m p - the 
mass of a proton, fi - molecular weight. The gas temperature 



1 http : //xmm. vilspa. esa. es/external/xmm_data_acc/xsa/ 
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9 = kT along the jet accounting for adiabatic and radiation 
cooling is determinated by the equation: 

3 d6 dn ( re 



= 0, 



(2) 



where e stands for specific energy losses by radiation per 
one particle. Introducing new variables r = 9/9o, x — r/ro, 
e — \inj one may write the equation in a dimensionless 
form: 

<h_ _ _4r _ J(t;6 ) 
dx 3 x x 2 ' 



(3) 



where a is a coefficient depending on physical parameters of 
the jet Mj and ro: 
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Qjm p v 2 r 9 



(4) 



The introduced value J(9) is the total emissivity of the 
plasma, and radiation losses integrated over all energies, 



(5) 



J{9) = / Je(9)de. 



The value J (8) is tabulated in the model of hot, optically 
thin plasma APE C/APED, which use s the atomic database 
ATOMDB vl.3.1 (ISmith et alJlioOlal . ISmith fc Brickhousd 
l200d , ISmith et al.ll2001bl ). ATOMDB provides an improved 
spectral modelling capability through additional emission 
lines, accurate wavelengths for strong X-ray transitions, and 
new density-dependent calculations (for more details see 
|http: / / cxc.harvard.e du / atomdb / ) . Relativistic boosting is 
taken into account in this model. The model provides a 
possibility of separate analysis of both line and continuum 
emissivities. 

To calculate the multi-temperature jet spectrum we di- 
vide the jet in 20 equal parts in the section ro — 4ro, the 
gas temperature at ~ 4ro drops belo w 100 eV, where th e 
thermal instabilities begin to operate (Ko tani et al.|[l996l ). 
One needs to consider many (at least 20) mono-temperatute 
parts of the jet, because the radiative cooling begins to dom- 
inate the adiabatic cooling at temperatures below ~ 2kev 
and has a strong peak at ~ 0.6 kev. The total spectrum of 
the jet is calculated using the formula: 



F e = ^NiJ e (e;) 



(6) 



where J e is the plasma emissivity and Ni is the spectrum 
normalisation of i-th part of the jet: 

1 



Ni = 



4tyD 2 



n e riHdV, 



(7) 
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wher e D is a distance to SS 433 (5.5 kpc, iBlundell fc Bowlerl 
2004). The integration is carried out over the jet portion 
volumes Vj. The average gas temperature of the i-th part of 
the jet 0* is estimated as follows: 

The main parameters which define the physical state of 
the jet are gas density and temperature #o(ro) at the visi- 
ble jet base. The density depends on the mass loss rate Mj 
and the jet (half) opening angle As we discussed in the 
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Figure 1. Measured ratio of the blue jet FeXXV and FeXXVI 
line fluxes (Table 1), derived with the simple local continuum fit, 
and the model ratio as depending on the gas temperature at the 
jet base. The highest temperature, $o ~ 12 keV, is obtained in 
orbit 707 (ip ~ 0.0), a slightly lower temperature in orbit 610 
(ip ~ 0.8); in orbit 609 we find the temperature 8o ~ 8 keV with 
notably bigger uncertainty. 



Introduction, the last two parameters were relatively well 
measured in SS433 in previous studies not only in X-ray 
range. The opening angle was particularly very well mea- 
sured in optical studies. To make our analysis more certain 
we do not try to determine the jet kinetic luminosity (the 
mass loss rate) from the XMM spectra, the more so, that one 
has serious problems with an interpretation of c ontinuum at 
~ 7 keV (Bri nkmann. Kotani. fc Kawai 2005|). Therefore, 
we adopt the parameters Lkin = 10 39 erg/s and 2-dj = 1.5° 
as already known. 

From all previous X-ray studies we may expect the gas 
temperature at the base is in the range of 0o(ro) =8 — 20 keV 
(see Introduction), and the visible jet base ro ~ 10 12 cm. 
The last estimate is quite substantial, it was argued above, 
that both the d onor star and the outer disc rim have about 
the s ame size (jCherepashchuk et al.l 120051 : iFilippova et al.l 
2006). The partial eclipses of the X-ray jets by the donor 
observed in SS 433 prove the estimate of the X-ray jet base 
ro- 

In the thermal jet model we adopt the solar abun- 
dance of elements, however, Nickel is 10 times overabun- 
dant (|Brinkmann. Kotani. fc Kawail l2005t ). Section [5] dis- 
cusses this in detail. 



4 JET PARAMETERS IN IRON LINE 
DIAGNOSTICS 

In this section we estimate the jet parameters 9o and 
ro using the fluxes and the flux ratio of the iron 
lines FeXXVKc and FeXXVILy a dKotani et all Il996l : 
iBrinkmann. Kotani. fc Kawail l2005h and adopting a sim- 
ple local PL continuum. These FeXXV and FeXXVI lines 
formed in the approaching jet are the strongest and well 
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Figure 2. Dependence of the derived gas temperature at the 
jet base on the variations of the visible base position tq (in the 
simple local continuum model) and mass loss rate Mj for orbit 
707 observation (ip = 0). The solid curve is for Mj = 3.3x 10 19 g/s 
(Lkin = 10 39 erg/s), the dashed grey curve is for 1.5 X 10 19 , the 
solid grey curve is for 5 X 10 19 g/s. Two horizontal dashed lines 
mark the temperature limits (Table 1). 




Figure 3. Determination of the blue jet visible base ro from the 
FeXXV line flux and the temperature derived (Table 1) for the 
observation in the precessional phase ip = 0. To show confidence 
levels for ro (two solid grey lines) we used triple errors of the 
line flux. Horizontal dashed line corresponds to the mass loss rate 
Mj = 3.3 x 10 19 g/s (L kin = 10 39 erg/s). 



resolved in the XMM spectra. The resulting line fluxes de- 
pend on the underlying continuum model. Therefore we use 
the simplest continuum to analyse here how stable are the 
derived parameters depending on Mj and on the errors of 
the measured line fluxes. We do not examine here different 
contin uum models, because iBrinkmann. Kotani. fe Kawail 
|2005l ) have studied these XMM spectra using many ver- 
sions of continuum and concluded that "the high sensitivity 
and wide bandpass of the XMM-Newton instruments rule 
out any of the simple continuum models used previously" . 
Instead in the section below we find clear evidences for a 
certain additional continuum component (reflection), study 
the whole spectrum with these two models, thermal jet and 
reflection, and correct the values of 9o and ro. 

In our model the iron line flux (as it follows from (HUE}) 
may be written as: 

M 2 

Fix J^-J-tp^Mj). (9) 

iljTo 

Here tp(ro,Mj) is a function depending on the gas cooling 
details and it does not strongly depend on ro and Mj. At 
the same time the line ratio FeXXV/FeXXVI depends on 
the gas temperature at the base 8o and its derivative d8/dr, 
which is determined by (|3]) . Modelling of the iron line fluxes 
gives us information on the temperature and visible jet base 
ro- 

To determine the line fluxes in the observed spectra 
we model the continuum with power law, and the lines as 
Gaussian blueshifted to z = —0.105 and z = —0.095 for 
precessional phases -0 = and ij) — 0.8 respectively. These 
blue shits were found in our fits, and they are in approxima- 
tive agreement with those expected from the precessional 
ephemeris. For better fits of these two lines with the con- 
tinuum we discarded the spectral parts below < 7 keV and 
8.0 — 9.5 keV, where there is a complex of emission lines 
Ni XXVII, Ni XXVIII, FeXXVIL/3 and FeXXVK/3. This 



spectral interval also includes a broad abso rption feature 
(an ab sorption edge), which was reported bv lKubota et al.1 
l|2007l ). In this simplest continuum model the line fluxes are 
well determined. Table 1 presents the fluxes of these two 
blueshifted lines and their ratios along with errors. 

In Fig.[T] we show the observed flux ratios 
Fe XXV / Fe XXVI along with the model ratio depend- 
ing on the gas temperature at the jet base. We check how 
stable is the temperature derived relating to the variations 
of Mj and ro. Fig. [2]shows the temperature as depending on 
ro at different mass loss rates Mj for the observation taken 
in the precessional phase ip = 0. These two parameters do 
not change the base temperature notably. We find the same 
for the observations in the precessional phase tp = 0.8. One 
may expect this result because the line flux © depends 
on Mj and ro in the same manner for both lines. Fig. [2] 
shows that to keep the same temperature at the visible jet 
base 0o (ro) (the same iron line ratio), with bigger Mj we 
have bigger ro to have about the same gas density no(ro). 
The function <p(ro,Mj) depends on the gas density, which 
is n oc Mj/r 2 . 

There is also some interplay between the estimates of 
#o(ro) and ro when we fix the mass loss rate in the jets. 
The emissivity of the FeXXV line peaks at 9 ~ 5 — 6 keV, 
and the total line flux is oc n 2 . When one finds 9q from the 
Fe XXV / Fe XXVI line ratio and determines the jet base ro 
from the FeXXV line flux, one obtains the gas density at 
the base no(ro). The higher the temperature 9o derived with 
regard to the 5-6 keV value, the lower the gas density in 
the emissivity peak jet region, the less ro we obtain to match 
the observed line flux. This effect is seen in Fig. [2] and it is 
not strong. 

Thus, we may find the temperature at the visible jet 
raduis #o(ro) from the Fe XXV /Fe XXVI flux ratios, and 
then we may find the base (ro) from the observed flux of the 
blue shifted Fe XXV line and the derived temperature. The 
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Table 1. Observed fluxes and derived blue jet preliminary parameters in the simplest continuum model. The columns present precessional 
ip and orbital tp phases of the observations, line fluxes in units 10 -4 ph/scm 2 , the flux ratios 1Z, the temperature at the jet visible base 
00 and the base ro. 



i> 


V 


F(FcXXV) 


F(FeXXVI) 


iz 


00, keV 


ro, 10 11 cm 


0.04 


0.28 


5.5 ±0.4 


1.95 ±0.3 


2.8 ±0.5 


12 ±2 


2.8 ±0.2 


0.84 


0.34 


4.9 ±0.7 


0.8 ±0.4 


6.1 ±3 


8± 4 


2.6 ±0.3 


0.84 


0.49 


5.9 ±0.2 


2.0 ±0.2 


3.0 ±0.3 


12 ± 1 


2.6 ±0.1 




Figure 4. A sketch of the funnel and the blue jet visibility at two 
precessional phases considered, ro (r' ) is the distance between the 
base of the visible jet and the top of the cone (the bottom of the 
funnel). r\ (rV) is the distance between the base of the visible 
funnel wall and the top of the cone. 



values of 0o(ro) and ro found with the simple continuum 
model are presented in Table 1. 

Fig-E] shows how the base ro is determined from the 
FeXXV line flux. We used here the triple errors (Table 1) 
of the He-like blue jet line flux. The blue jet visible base we 
found is ro ~ 3 • 10 11 cm. The equation ((9]) provides us with 
useful and obvious scaling formula, ro oc Mf /fij, and Fig. [3] 
confirms the scaling. With different jet parameters satisfying 
this scaling and with about the same gas temperature at 
the visible jet base, we obtain about the same thermal jet 
spectra. 

Fig.[l]shows a sketch of the jet and the supercritical ac- 
cretion wind (funnel) geometry at two precessional phases 
which we consider, ip « (orbit 707) and ip ~ 0.8 (orbits 
609, 610). They are marked in the figure as ro and r' re- 
spectively. We adop ted the funnel opening angle $/ = 25° 
(jOhsuea et al.f2 005) in the figure. Our estimates of r and r' 
are the same within the errors (Table 1), however the result 
depends on the fluxes of the Fe XXV line, i. e. on the natural 
variability of SS 433 jet activity. In the second observation 
at the ip ~ 0.8 the iron line is brighter. 

The estimates of 6o(ro) and ro obtained in this section 
were derived using the simplest continuum model. In this 
model we discarded the spectral part below < 7 keV and 
the region of 8.0 — 9.5 keV, w hich includes a broad absorp- 
tion edge l|Kubota et al.| [2007). This absorption feature may 
distort the Fe XXV / Fe XXVI line ratio making it higher, as 
a part of the Fe XXVI line profile may be insi de of the broad 
edge. The edge profile is not established yet l|Kubota et al.l 
120071 ). and we can not rule out that real temperature 80 (ro) 
is somewhere higher than the one we derived here. Our es- 
timate of ro is not affected by the absorption edge, because 
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Figure 5. Observed spectrum of SS433 (orbit 707, ip ~ 0) with 
the thermal jet model components and the additional reflection 
model in two spectral regions 5.0—7.0 (top) and 7.0-10.0 keV (bot- 
tom). The total model spectrum is shown with solid line. The re- 
flection model explains well the fluorescence iron line at 6.4 keV, 
the recombination iron K a line of Fe XXV at 6.7 keV (top) and 
the iron absorption edge llKubota et alj|2007ri at ~ 8 - 9 keV. 



it is based on the FeXXV line. We use these estimates of 
@o(ro) and ro (Table 1) as a first approximation. Below we 
find more accurate values of these parameters using a more 
complex model of the continuum, which describes the whole 
spectrum of SS 433. 



5 ADDITIONAL COMPONENTS IN SS 433 
SPECTRA. REFLECTION MODEL 

As it was noted in the Introduction, there is a problem of 
four orders of magnitude, missing in X-ray luminosity of 
SS 433. One may expect to find the original luminosity of the 
object L x ~ 10 40 erg/s, instead, we observe L x ~ 10 36 erg/s. 
This may mean that direct radiation is blocked inside the 
hydrodynamical funnel, where the jets are accelerated and 
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collimated l|Fabrikall20d3) . Taking into account such an ex- 
treme original luminosity we may expect to observe a reflec- 
tion of a part of the radiation at the outer funnel walls (the 
wind). In Fig. [4] we show the probable reflecting area on the 
wall at r > r\. 

There are strong observational indications that point to 
an existence of a reflection component in the X-ray spectrum 
of SS433. Firstly, the presence of a strong fluorescent line 
of quasi-neutral iron (6.4 keV), which can be a result of 
the reflection. Secondly, in the spectral region ~ 8 — 9 keV 
there is a broad abs orption feature (a trough), detected by 
iKubota et aD (|2007l ) . The trough may be an iron absorption 
edge arising in reflection, its position corresponds not to the 
neutral pl asma, but the some intermediate stages of iron 
ionization (Kubota ct al. 2007). This leads to an assumption 
that the reflecting matter in partly ionized. That is why we 
decide to fit the observed spectrum of SS 433 with the multi- 
temperature jet model plus reflection in the partly ionized 
material. 

For the reflection component we use the REFLIOr4j 
mode l of X-ray ionized re flection (jRoss. Fabian, fe Yound 
1 19991 ; iRoss fe Fabianl [20051 ). The illuminating radiation is 
assumed to have a cut-off power-law spectrum Fe = 
AE~ F+1 exp(— E/E c ), where the cut-off energy is fixed at 
E c = 300 keV. It is assumed that the surface is irradiated 
with the X-rays so intense that its ionization state is deter- 
mined by the ionization parameter, 

£ = (10) 

n 

where F to t = / FecIE - the total illuminating flux, n - the 
surface gas density. The amplitude A in the total illuminat- 
ing flux is chosen to provide a desired value of the ionization 
parameter £. In the model the irradiated gas is Compton 
thick and has a constant density. The illuminating radiation 
flux has a power-law spectrum in the model, and the photon 
index F is a model parameter. 

Fig-El shows the observed spectrum of SS433 with 
the thermal jet model and additional reflection model. It 
presents two regions 5.0-7.0 keV and 7.0-10.0 keV fit- 
ted with the same model parameters. Both the reflection 
emission feature at 6.4 keV and the iron absorption edge 
at ~ 8 — 9 keV are reproduced well in the reflection 
model. Fig. discovers one more an evidence of the reflec- 
tion in SS433 spectrum, it is an emission feature of iron 
FeXX VKq at 6.7 keV at zero velocity. The REFLION 
model (|Ross fe Fabianl [2005) considers the Ka recombina- 
tion lines of FeXXV-XXVI (near 6.7 and 7.0 keV) and the 
fluorescence lines of Fe VI-XVI (near 6.4 keV). The fluores- 
cence features never been explained before in SS 433 spec- 
trum, this convinces us to use the reflection model as an 
additional (to the thermal jet) component. In the reflection 
model fit we estimate the ionization parameter at the re- 
flecting surface of £ ~ 300 and the illuminating radiation 
spectral index of T ~ 2. The very good continuum fit allows 
us to estimate more accurately the iron line fluxes in the 
blue jet and respectively the visible jet base temperature, 
it is #o(ro) ~ 17 keV. Below in this section we study these 
parameters in more details using the whole spectrum. 

1 http : //heasarc . gsfc .nasa.gov/docs/xanadu/xspec/models/ 
reflion.html 




Figure 6. Observed spectrum of SS433 (orbit 707, if) ~ 0) 
with two thermal jet models. The first model (8o(ro) = 17 keV, 



Njj = 1.5 X 1CH 



) is shown with solid thin line, individual 



parts of the multi-temperature jet are shown with dashed lines to- 
gether with the model's residuals in bottom. The model produces 
well the blue jet lines in the region 7—9 keV. The second model 
(00 (ro) = 14 keV, N H = 1.0 X 10 22 cm" 2 ) is shown with grey 
thick line, this model is intended to fit the lower energy range, 
however it fails to produce observed line fluxes of the blue jet. 
The strong excess of the observed continuum radiation in the re- 
gion > 2.5 keV may be fited well with the reflection additional 
component (section 5). 



The harder part of the spectrum > 5 keV is not sensi- 
tive to the absorption value Nh, from the other hand, the 
reflection component we found weakens notably in softer en- 
ergies, < 4 keV. That is why to study the absorption value 
Nh, in a next step we analyse the whole spectrum of SS 433 
with the thermal jet component only. 

In Fig.|5]we show the observed spectrum of SS 443 along 
with the thermal blue jet models. We only used in our figure 
the blue jet, because adding the red jet can not in effect 
solve the continuum problems both at higher energy an d 
even at lower energi es (|Brinkmann. Kotani. fe Kawaill2005f ) . 
iKotani et al.l (|l996t ) have found that the receding red jet is 
not visible in low energy ranges (in the ranges < 4 keV) in 
the precessional phases ip ~ 0, when the accretion disc is 
the most open to the observer. They interpreted this fact 
via an existence of absorbing gas in the precessing plane 
perpendicular to the jets (the accretion disc plane), which 
only blocks the outer cooler parts of the red jet. The hotter 
parts of the red jet located closer to the central machine, 
are visible, the extended material does not obscure them. 
The extended pre cessing disc-li ke outflow in SS 433 has been 
discussed before (|Fabrikalll99"3h . 

There are two jet models in Fig. [6] in both mod- 
els we have adopted our standard jet parameters, Mj = 
3.3 • 10 19 g/s and 2$, = 1.5°. There are strong dis- 
crepancies between the observed and the model spectra. 
Whereas the jet lines may be well fitted with one of the 
model (the FeXXV and FeXXVI lines of the approach- 
ing jet), t he continuum divergence is dramati c. This was 
noticed bv lBrinkmann. Kotani. fe Kawail l|2005h . They con- 
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eluded that when the line spectrum is well reproduced 
by the thermal jet model, but the continuum is not. We 
tried to use our jet model w it h the parameters found by 
iBrinkmann. Kotani. fc Kawail (|2005l ). and confirmed their 
results. When the hard energy region of the continuum can 
be fitted well by adding the reflection component, the soft 
energy range is not. 

Producing the fits shown in Fig. [6] we tried to match the 
spectrum best possibly to illustrate the thermal jet model 
limitations. The first model fits well the iron line fluxes, 
however we observe a strong soft deficit. The neutral hy- 
drogen column density Nh in the fit is 1.5 x 10 22 cm" 2 and 
the visible jet base temperature is 9o(ro) = 17 keV in the 
model. Note this column density is in very good agreement 
with the interstellar e xtinction in S S433 found in optical 
and UV observations l|Fabrikall2003 ). To diminish the soft 
deficit one must take less temperature at the jet base. The 
best fit we can find in this way is shown in Fig. HJ] as the sec- 
ond model (N H = 1.0 x 10 22 cm -2 , 6 (r ) = 14 keV). It has 
not such a big problem in the low energy region, however 
the blue jet iron line fluxes can not be reproduced in this 
model, they are at least ~ 1.5 times less than the observed 
fluxes. Besides that, the softest energy spectrum (< 1 keV) 
is not descibed well even in the second model. We conclude 
that the thermal jet model can not describe simultaneously 
the line fluxes and the soft continuum. 

We have found that the thermal jet model produces 
about the same spectra when we keep the same value of 
M? /Q,jTo = const. This scaling formula is valid in the first 
approximation only, ip(ro,Mj) in ([9]) does not strongly de- 
pend on ro and Mj (see previous section). For instance, when 
we compare the thermal jet model shown in Fig. [6] by thin 
solid line (Mj = 3.3 • 10 19 g/s and 2dj = 1.5°), with another 
model, where we take these two parameters twice as bigger 
(in a limit of small angles Qj oc i? 2 ), we obtain the same 
model spectrum when we increase the value of Nh by 17%. 
Both the iron line fluxes and their equivalenth widths are the 
same in these two models. However the line profiles appear 
too broad in the second case, at least during this (orbit 707) 
observation the jets were notably narrower than 2i?j = 3.0°. 

It is very importan t to note here that 
(Brin kmann. Kotani. fc Kawail [2005) have concluded 
that the thermal jet model in any combination can not pro- 
duce the observed continuum. We confirm their conclusion 
here, and we may state the problem more concrete. When 
one reproduces well the continuum around the blue jet 
iron lines (~ 6 — 10 keV), one can measure the line fluxes 
without strong uncertainties. One can determine then the 
jet visible base temperatute and estimate the value of 
using the thermal jet model. We have found 
in this way that the thermal jet model can not produce the 
lower energy continuum in SS433. However the strongest 
divergence is observed in higher energy region, where about 
50 % of the flux is formed not in the thermal jet (Fig.0. We 
have discussed above the clear evidences of the reflection in 
this spectral region. 

We know nothing about the incident flux and its spec- 
trum in SS433, so we divide the whole spectrum into four 
parts (0.8-2.0, 2.0-4.0, 4.0-7.0 and 7.0-12.0 keV) and find 
the reflection model parameters independently in each part, 
but with the same jet parameters and the same IS absorp- 
tion Nh- 



The red jet spectrum has to notably contribute to 
the low energy reg ion, however, it is not observed there 
(|Kotani et alj fl996). but it is observed at higher energies 
(red iron lines). The lack of the red jet at low energies 
(< 4 keV) may be bound not only to its specific visibility, 
but with an underestimated extinction in the whole spec- 
trum of SS 433 as well. From the XMM spectra available we 
can not realise clearly which effect is correct, an obscuration 
of a cooler portion of the red jet or an underestimated extinc- 
tion in SS 433. Therefore we produce two fits, first without 
the red jet cooler portion in the range 0.8 - 4.0 keV (case A), 
and second with the red jet in the whole spectrum with in- 
creased Nh (case B). We find below that Nh values derived 
in both cases are about the same and this does not change 
our conclusions on the additional source of X-ray radiation 
in SS 433. 

We estimated the column density Nh of the neutral IS 
absorbing material in the range 0.8-2.0 keV. In this range 
a strong excess of radiation is observed (Fig©, which does 
not practically depend on the thermal jet model parameters. 
Hereat we added a BB component, which might be an indi- 
cation of the supercritical accretion disc funnel (a wind) ra- 
diation (see Section|5J. In Fig. [7] we show the observed spec- 
trum of SS433 along with the model fits 'jet-plus-reflection' 
in the four spectral ranges. Table[2] presents the model pa- 
rameters. The thermal jet model is identical in all ranges, 
yet the reflected component was fitted independently in each 
range (however there is no reflection in low energies). 

The blue jet parameters were found in the high en- 
ergy range 7.0-12.0 keV matching the model fluxes and 
flux ratios of the He- and H-like iron lines to the observed 
values. We do not see any need to use non-solar abun- 
dances of the elements, but Ni is clearly overabundant. 
IBrinkmann. Kotani. fc Kawail l|2005l ) found the Ni abun- 
dance to be about 8 times solar. Our fits show that Ni has 
to be overabundant by a factor 10, since the Ni line is inside 
the broad absorption feature (Fig.[7]d). The red jet param- 
eters were found in the range 4.0-7.0 keV, where the red 
iron lines are seen clearly (Fig.[7]c). 

The jet parameters presented in Table[2] are more ac- 
curate than those presented in Table 1, where we used the 
simplest continuum model. In the precessional phase ip ~ 0.0 
we find that the blue jet visible base is of ro,& ~ 2 ■ 10 11 cm 
and the gas temperature in this place is do ~ 17 keV. The 
red jet base is of ro, r ~ 7 • 10 11 cm, which is naturally ex- 
plained by the obscuration of the red jet by the accretion 
disc (the wind) body in this orientation of the disc. In the 
precessional phase ip ~ 0.8, where one may expect an ob- 
scuration of a bigger portion of the blue jet we find a lower 
temperature, 9o ~ 15 keV. However, ro,6 is about the same 
in both precessional phases. The spectra taken at ip ~ 0.8 
are notably noisier than the one taken at ip ~ 0.0. 

Finding accurate jet parameters is not the main pur- 
pose of the present paper. The results depend on the natural 
variability in the jet activity, besides that the physical pic- 
ture may be more complex (possible gas clumping, heating 
effects due to kinetic energy dissipation or the funnel radia- 
tion). However, the derived jet visible base ro,b ~ 2 ■ 10 cm 
is too small, while it has to be about 10 12 cm to satisfy 
the observed eclipses by the donor star. Th e inner hot- 
ter portions of the jet are partially eclipsed (|Kawai et al.l 
ll989l : lFilippova et alj|2006h . The jet's length is scaled © as 
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Figure 7. SS 433 spectra (orbit 707, if) ~ 0) with the model 'jets-plus-reflection' fit in four spectral ranges for case A, with zero red jet 
contribution in the ranges 0.8-4.0 keV (Tabled- The brightest jet components of the multi-temperature model are shown by dotted 
lines for the blue jet and dash-dotted lines for the red jet, the total model spectrum is shown by solid line, a) 0.8-2.0 keV range, where 
we added a black body spectrum to describe a soft excess observed at ~ 1 keV (dashed line), b) 2.0-4.0 keV, where the reflected 
spectrum (dashed line) contribution weakens notably at lower energies. It is indicated by some excess in the observed spectrum at 
~ 2.2 — 2.6 keV. c) 4.0 — 7.0 keV, where the reflection well describes both the continuum radiation, the fluorescent line of quasi-neutral 
iron and the recombination FeXXV line near 6.7 keV. Some excess at lower energies of the observed continuum is related to the change 
of the incident radiation spectral index (idem as in the previous range), d) 7—12 keV range is very well described by the reflection both 
in the continuum and in the broad absorption feature dKubota et al.ll2007Tl . 



ro tx Mj/Zlj, and one may scale the length to about the 
expected value of ~ 10 12 cm in the thermal jet model which 
we use. 

The reflection model produces a very good fit in the 
high energy range 7.0-12.0 keV (\ 2 = 1.05), it describes 
well both the whol e continuum and th e broad absorption 
feature detected bv lKubota et all (|2007l ). The feature is lo- 
cated between the blue shifted FeXXVI and Ni XXVII lines 
and partly covers these lines (Fig.[7]d). The parameters of 
the reflection are found, the ionization parameter of the ir- 
radiated surface is £ ~ 300 and the illuminating radiation 
index is V ~ 2. This means a flat (EFe) spectrum of the 
incident radiation in this range. 

In the range 4.0-7.0 keV we have an equally good 
agreement between the observed spectrum and the model. 
The reflection features at 6.2-6.8 keV are well produced. 



However, in this range there is a small continuum divergence 
at ~ 4 keV, probably due to the fact that the illuminat- 
ing radiation index changes with energy. We find the index 
r ~ 1.6 in the range, the ionization parameter is about the 
same as at higher energies, £ ~ 240. 

In the range 2.0-4.0 keV we observe the same effect 
is stronger. The model has a some deficit at 4 keV and an 
excess at 2 keV. An inclusion of the reflected component is 
essential in the range, otherwise one can not produce the ob- 
served continuum at 3-4 keV. The incident radiation has to 
have an index F < 1 in the range. This is confirmed by a too 
strong emission feature at w 2.5 keV in the reflected model. 
Similarly to the previous spectral range we may expect that 
the incident radiation index is changed with energy, and here 
a steeper incident spectrum is expected. The model we use 
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Table 2. Results of the spectral fits in our four spectral ranges. The upper part corresponds to the 
observation of orbit 707 (ip ~ 0.0), the bottom part, - to the orbit 609 (ip ~ 0.8). Case A is a model with a 
zero red jet contribution in softer (0.8 — 4.0 keV) ranges, case B with the red jet included. The jet parameters 
are Bq ss 17 keV, r b Kb 2 ■ 10 11 cm, ro, r ~ 7 ■ 10 11 cm in the precessional phase ip ~ 0.0, and 6o ~ 15 keV 
and the same visible jet base in ip ~ 0.8. The lines present the ionization parameter (1101 (£, erg/cm s), the 
photon index of a power law illuminating radiation T, the BB temperature in the soft X-ray range 8 bb (keV), 
the normalisation of the BB spectrum Normj,;, = L^g/D^Q (L39 is luminosity in units 10 39 erg/s, Diq is a 
distance in units lOkpc), and column ISM density Nh in units 10 22 cm -2 . Nh is equal in each range as 
found at lower energies. 



Range, keV 




7-12 


4-7 


2-4 


0.8-2 




A 


300 ±5 


240 ± 30 


si 300 


- 




B 


300 ±4 


250 ± 30 


m 300 


- 


r 


A 


2.00 ± 0.08 


1.56 ±0.05 


< 1 






B 
A 


2.00 ±0.09 


1.57 ±0.04 


< 1 


0.107 ± 0.001 




B 








0.100 ±0.001 


Norm^f, 


A 








0.032 ±0.002 


B 








0.125 ±0.005 


N H 


A 








1.54 ±0.01 


B 








1.85 ±0.01 


X 2 (d.o.f) 


A 
B 


1.05 (441) 
1.13 (441) 


1.28 (535) 
1.13 (555) 


1.39 (396) 
1.55 (397) 


1.39 (242) 
1.58 (242) 




A 


340 ± 80 


300 ± 90 


300 




B 


310 ± 10 


300 ± 90 


« 300 






A 


2.2 ±0.2 


1.56 ±0.08 


< 1 




r 


B 

A 


2.0 ±0.2 


1.57 ±0.08 


< 1 


0.093 ± 0.004 




B 








0.097 ±0.003 


Norm;,;, 


A 








0.026 ± 0.003 


B 








0.068 ± 0.005 


N H 


A 








1.28 ±0.01 


B 








1.58 ±0.01 


Xl (d.o.f) 


A 
B 


1.36 (189) 
1.39 (190) 


1.1 (332) 
1.08 (322) 


1.01 (388) 
1.05 (388) 


1.55 (234) 
1.05 (234) 



is restricted by the value F = 1.0. We find then at F = 1.0 
the ionization parameter £ ~ 300 (with a notable scatter). 

In the soft range 0.8-2.0 keV the reflected component 
is not detected. Even if there is a reflection in the soft range, 
it has other parameters than those found at higher energies. 
From the 7.0-12.0 keV band to the 2.0-4.0 keV band, we 
found a gradual decrease of the spectral index of the illu- 
minating radiation, - from F « 2.0 to F < 1 at about the 
same ionization parameter, £ ~ 300. However, an additional 
component related to the thermal jet one, is indeed present 
in the soft energy range (Fig. [6] and Fig.[7]a). Choosing the 
simplest way we add a BB component to describe the soft 
excess detected in the range. 

The temperature of the BB component is 6 bb ~ 0.11 keV 
and it is identical in both our cases, A (without the red jet) 
and B (with the red jet). The BB component dominates at 
energies < 1 keV and it is true again in both cases. When we 
do not include the red jet in the fit, the BB luminosity is ~ 4 
times smaller compared to the one with the red jet. However, 
this is rather a result of the increased Nh in the latter case 
(B). The soft ~ 1 keV radiation is very sensitive to the Nh 
value, therefore we can not estimate the real luminosity in 
the soft excess from the XMM spectra. 

In both cases we find the IS extinction value to be about 
the same, it is (Table[2j N H « 1.54 x 10 22 cm -2 in case 
A and N H ~ 1.85 x 10 22 cm -2 in case B. The interstellar 



extinction in SS 43 3, found in optical and UV observations 
ilDolan et alj 119971), is A y = 8.0 - 8.4. This corresponds 
jpredehl fc Schmitt|[l995l ) to N H = (1.4 - 1.5) x 10 22 cm -2 , 
which is in good agreement with the extinction found in 
X-rays. 

Spectral residuals in the soft energy range are bigger 
than the ones we have in other ranges. The most prominent 
negative feature which is not produced in the model fit is 
at 1.2 — 1.3 keV (Fig. 0a). However, such a feature was pre- 
dicted to be present in the spectra of ULXs, if they possess 
a supercritical accretion di sc similar to the one in SS433, 
but ob served nearly face-on. iFabrika. Abolmasov. fc Karpovl 
l|2007| ) (see their Fig. 2) have simulated spectra of the multi- 
colour funnel model (MCF), where they introduced an Lc 
edge of O VIII, blueshifted to the SS 433 jet velocity. The 
absorption features arise inside the funnel. The model spec- 
tra being reduced with standard XSPEC models with IS 
absorption (Nh ~ 10 22 cm -2 ) produce fake broad emis- 
sion/absorption features. The strongest one is an absorption 
feature at 1.2-1.3 keV. 

Our model fits to the spectra taken at the precessional 
phase ip ~ 0.8 give us approximately the same results (Ta- 
ble^ . We find a little lower temperature and less luminosity 
of the soft BB component. In the sketch presented in Fig. 2] 
the probable reflecting region is shown at r > n at ip ~ 
and r > r[ at ip ~ 0.8 on the visible outer wall of the funnel. 
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Figure 8. Additional components in X-ray spectrum of SS 433 
(the thermal jet component is not shown). The illuminating in- 
cident radiation found in reflection model is shown by solid 
lines (arbitrary scaled to ELe = 10 40 erg/s) in three energy 
ranges: 2.0-4.0 keV (r = 1.0), 4.0-7.0 keV (r = 1.6) and 7.0- 
12.0 keV (r = 2.0). The soft BB component (9 bb = 0.1 keV, 
L x ~ 3- 10 37 erg/s) detected in the 0.8-2.0 keV energy range is 
shown by solid line. There are two originally the same flat spectra 
after absorption in ABSORI model with N H = 10.0 X 10 22 cm -2 , 
£ = 150 (grey dash-dotted line) and N H = 6.0 X 10 22 cm" 2 , £ = 3 
(bold dash-dotted line). There are two MCF (multi-colour funnel 
model) non-absorbed spectra with ri n = r\ = 1.74- 10 11 cm (bold 
dotted line) and = 5 ■ 10 10 cm (grey dotted line). 



If the soft excess is the visible funnel wall's own radiation 
(plus a possible reflection of a soft funnel radiation), one 
may expect less luminosity of the wall due to the visibility 
conditions in the precessional phase ip ~ 0.8. 



6 DISCUSSION 

In this section we first discuss the structure of the super- 
critical accretion disc wind based on the wind approach by 
IShakura fc Svunvaevl (| 19731 ). and estimate the wind and the 
jet funnel parameters for SS 433. Next we discuss the re- 
flected spectrum's origin and possible interpretations of the 
soft excess. 



6.1 The expected parameters of the funnel in 
SS433 

The limiting Eddington luminosity is Le — AttGMc/k = 
1.5 • 10 39 mio erg/s, corresponding limiting mass accretion 
rate is Me = 48nGM/ck = 2 ■ 10 19 mio g/s, where the 
Thomson opacity for a gas with solar abundance is n = 
0.35 cm 2 /g, and the black hole mass mio is measured in 
10 solar mass units. With an accretion rate at the outer 
disc boundary M a » Me, the disc becomes supercriti- 
cal at spherisation radius r sp . In its outer parts r > r sp 
the disc is subcritical, M(r) = M a , below the spherisa- 
tion radius r < r 3p the accretion regime is locally critical, 
M(r) ~ (r/r 3p )M a and there appears a strong wind. The 
total mass loss in the wind below the spherisation radius is 



M w (r) ~ M a (l — r/r sp ). The gas is blown away by the ra- 
diation pressure and the accretion rate is decreased linearly 
with radius. The outflow velocity is close to that of Kep- 
lerian, v w (r) ~ (GM/r) 1 ^ 2 , for the bulk of the wind it is 

Nl/2 



(GM/r sp ) 1/2 . Following IShakura fc Svunvaevl 
the spherisation radius is (also see lPoutanen et aLn f 



3kM a 
8ttc 



(11) 



This simple picture first drawn bv lShakura fc Svunvaevl 

l| 19731 ) neverthel ess explains th e main features that we ob- 
serve m SS433 (|Fabrikall200i ). It is also confirmed in the 
main description in radiation hydrodynamic simulations 
(|Ohsuga et al.l 12005). The simple model does not account 
for advection effects in the disc regions close to a black 
hole, however, this does not change the main picture. The 
M w ~ M a stays the same if M a » Me, even the emerg- 
ing spectrum of the supercritical dis c does not change no- 
tably when advection is included rtPoutanen et ail 120071 : 



fl 

iTakeuchi. Mineshige. fc Ohsugdl2009T T 

Adoptin g the mass exchange rate in SS 433 M w ~ M a ~ 
10~ 4 Mq /y (jShklovskiill 198lt Ivan den Heuvelll98ll ; iFabrikal 
2004) we find the spherisation radius r. 
and the wind velocity v w (r sp ) ^ 
where m_4 is measured in 10 -4 Mq/j units. The wi nd ve- 
locity is directly measured in SS433 (Fabrika 2004), it is 
~ 100 km/s at the polar angle 90° (disc edge-on) and in- 
creases sharply to ~ 1500 km/s at 60° when the disc is the 
most open to observer. 

Inside the spherisation radius the disc luminosity is lo- 
call y Eddington, which g i ves t he total accretion luminos- 
ity (|Shakura fc Svunvaevl Il973h for SS433 L x ~ L E (1 + 
In M a / Me) ~ 1.0- 10 40 m 10 erg/s. Outside the funnel the ra- 
diation is thermalised (Lboi ~ L x ), the wind's photospheric 
radius and its temperature in the photosphere are estimated 
as 



10 10 m_ 4 cm 



4000(mi /m_ 4 ) 1/2 km/s, 



T p h = kM w /An cos $fV w ~ 
T ph = (L bo i/47vacos-dfrl h ) 1/4 



1 • 10 12 cm, 



. ■ 10 4 K. 



(12) 



Both the size of the hot body and its temperature are 
quite close to those we observe in SS433. Inside the fun- 
nel where the jets are formed, the photospheric radius (the 
funnel bottom) is r p h,j = nMj/QfVj ~ 2.5 • 10 9 cm, where 
Qf = 27t(l — cos 1?/) is the funnel's solid angle and #/ = 25°. 
The SS433 jet funnel is not transparent down to the black 
hole. Below r p h,j the funnel walls are probably not seen di- 
rectly by a face-on observer. 

A multicolour funnel model (MCF ) has 
been proposed for SS 43 3 dFabrika et all 120061 ; 
iFabrika. Abolmasov. fc Karpovl |2007| ) to estimate the 
emerging X-ray spectrum. They considered the opaque 
disc wind and estimated the gas cooling and its temper- 
ature inside the wind. In the inner parts of the wind the 
radiation pressure dominates and T(r) cx r -1 / 2 , in the 
outer parts the gas pressure dominat es and T(r) oc r~ . 
The supercritical disc simulations (|Ohsuga et all I2005T ) 
show that the T(r) oc r~ 1//2 is valid from r ~ 20 r 9 (below 
which the advective energy transport is important) to 
r ~ 200 r g ~ 10 9 cm, approximately the outer boundary of 
the computational domain. 

One may estimate the radiation trapping radius in 
the wind outflow comparing the dynamical time-scale 
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tmov ~ r/v(r) and the photon diffusion time-scale 
t esc ~ n(r)r 2 aT/c, where n(r) = M w /4tt cos i}frnpr 2 v(r), 
ot is Thomson cross-section, and mp is proton mass. 
The critical trapping radius in the wind (where t esc = 



IS Ttr 



(2/cos# f 



2 r s „ . One may ex- 



imov ) 

pect (|Fabrika. Abolmasov. fc Karpovll2007l ; iPoutanen et al.l 
I2007I ) that radiation is trapped within the wind up to 
r ~ 2r sp . In this region the radiation transfer is local. In 
outer regions r » 2r sp the radiation transfer inside the 
wind is global and radiation cools the wind (T(r) oc r _1 ). 

This simple approach allows us to estimate (see be- 
low) the soft excess in the spectrum observed in the low 
energy range (Fig. 0a) . Adopting the wind temperature of 
T p h ~ 60000 K at the outer photosphere one may expect 
(T(r) oc r~ ) the temperature of T ~ 0.5 keV at the level of 
spherisation radius. 

Photons entering the funnel from the funnel's walls have 
multiple scatterings on the walls before they may finally 
escape the funnel. The probability to escape from the infinite 
funnel after one scattering is P = 1 — cos $ / . The probability 
for a photon to remain in the funnel after N scatterings 
is (1 — P) . If we consider that 90% of photons escaped 
the funnel (cos N tf f =0.1), we find the average number of 
scatterings AT(90%) ~ 23, when ti s = 25°. The number 
of scatterings in SS433 funnel (r p h/r p h,j ~ 1000) is not 
significantly smaller than that in the infinite funnel. 

At the wall temperature of less than a few keV, as 
we may expect, the absorption is more probable than elec- 
tron scattering. This means that the photons are re-radiated 
back to the funnel with the local wall temperature T(r). 
The multicolour funnel has to have an integral spectrum 
l|Fabrika et al.l 120061 ; iFabrika. Abolmasov. fc Karpovll2007h 
with numerous edges, where the strongest edges are from the 
most abundant elements like iron and oxygen. The underly- 
ing spectrum with the edges formed at about zero velocity 
is absolutely necessary to produce the observed stable jet 
outflow velocity Vj 0.26c in the line-locking mechanism 



i|Shapiro, Milgrom. fc Reeslll986l ; iFabrikal 12004 ) 



6.2 The reflected and black body component 

Both the reflected component which we detect in three en- 
ergy ranges (Fig. [7] Table [5J), and the BB component in the 
soft energy range 0.8-2.0 keV, may be formed at the outer 
visible part of the wall shown in Fig. [3] at r > r\. At the pre- 
cessional phase ip « we find ri ~ 0.87 ro « 1.74 ■ 10 11 cm 
(Table [2]| . The ionization parameter (|10[) found at the re- 
flecting surface is £ ~ 300 and it is approximately the same 
in all the three energy ranges. As we discussed above, the 
multiple scattering in the funnel is expected, so we may con- 
sider this value of £ as effective or averaged over all putative 
scatterings in the funnel. A study of the additional spectral 
components is out of the scope of this paper, we only make 
rough geometrical estimates here. 

We estimate an angle j3 at which the reflecting medium 
is exposed to the incident radiation to produce the derived 
ionization parameter (110p . (3 ~ £nr 2 /L x . Using n(r)r 2 = 
M w /4n cos $fmpv m , we find f3 ~ ^M w /Ati cos , &jmpv w L x ~ 
0.04 m_4 ?j^ q I/Jq 1 for £ = 300, where the wind velocity is 
measured in units of 2000 km/s and the funnel luminosity 
in 10 40 erg/s units. This corresponds to the funnel depth 
f3rx ~ 7 • 10 9 cm, which is close to the spherisation radius. 



The reflected component luminosity found in the energy 
range of 7- 12 keV is L re fi ~ 3.3 • 10 35 k erg/s (at a distance 
of 5.5 kpc), where the k factor is ratio of the total cone 
surface at r > n (Fig.[4j to the visible cone surface. At 
the precessional phase ^i*0, the ratio is k « 3.2, and the 
luminosity at r > ri is L re /; ~ 1.0 ■ 10 36 erg/s. In the range 
of 7- 12 keV the reflection albedo is mainly geometrical and 
we might expect the angle j3 ~ F re fi/Fi nc , where Fi nc is the 
incident flux in this energy range, if the reflecting surface 
is directly exposed to the hard radiation from its origin. 
Assuming that the illuminating spectrum is flat, with total 
luminosity of ~ 10 40 erg/s in the spectral range of 0.1 - 
20 keV, we find that its luminosity in the 7-12 keV range 
is F inc ~ 10 39 erg/s and the angle [3 ~ F ref i/F inc ~ 0.001. 
This gives us an estimate of the funnel depth, where the 
illuminating radiation is formed, /3ri ~ 2 • 10 8 cm. 

The funnel photosphere radius, r p h,j ~ 2.5 • 10 9 cm, de- 
rived from the observed values (Mj,Vj), is notably greater 
than the last estimate of the radius where the hard illu- 
minating radiation is formed. Indeed, the reflecting surface 
may not directly see the deepest funnel parts. There have 
to be multiple scatterings in the funnel where the harder 
(> 7 keV) radiation may survive absorption. 

In Fig. [8] we show the illuminating spectrum found 
above in the 4-12 keV energy range, arbitrarily scaled to 
the luminosity of 10 40 erg/s. We show the spectral slope in 
the 2-4 keV range T — 1.0, but it was noted above that it 
has to be steeper, F < 1. The figure also shows two originally 
flat absorbed spectra. These two absorbed spectra illustrate 
two limiting cases between which one may expect to find 
true picture. Even very thick (Nh ~ 10 23 cm -2 ) medium 
of higher ionization (£ = 150) can not depress the radiation 
spectrum below 3-4 keV. Lower ionization medium (£ = 3) 
does absorb the radiation effectively. 

Therefore we expect an existence of an absorbing (re- 
flecting) medium inside the funnel, which produces such 
a curved spectrum of radiation, illuminating the reflected 
surface. Note that having He- and H-like absorption edges 
at about zero velocity is a mandatory property of the 
funnel spectrum to produce the observed stable jet ve- 
locity (v-j x 0.26 c) due to the line-locking mechanism 
l|Shapiro. Milgrom. fc Reeslll986l : iFabrikailiool ). 

Inspecting the results of t he sup ercritical accretion disc 
simulations by lOhsuga et all (|2005T 1 we find that the walls 
at the level of ~ 10 8 — 10 9 cm consist of a low velocity gas 
with the temperature of T ~ 2 ■ 10 6 K and the density of 
n ~ (1 — 3) • 10 18 cm -3 . The walls are the source of strong 
wind, they supply the funnel with the gas. The funnel walls' 
photosphere with the temperature from T ~ 10 7 K at r ~ 
10 s cm to T ~ 10 6 Katr~ 10 9 cm, may be the source of 
the edged radiation spectrum. 

The nearly flat spectrum of the funnel (Fig.[7][8]), which 
is expected to illuminate the reflecting surface, extends at 
least from 4 to 12 keV. The supercritical accretion disc spec- 
trum is expected to be flat ()Poutanen et al.l [2007). however, 
it is supposed to go up to a few keV only. The Comptoniza- 
tion may extend and flatten the spectrum to higher energies. 
At an accretion rate of ~ 10 % of the Eddington rate, the 
local emitted spectrum in the inner region of the disc is 
approximated b y a diluted BB with the sp ectral hardening 
factor / ~ 1.7 (Shi mura fc Takaharal[l9 95). For the accre- 
tion rates close to the Eddington limit, / increases and may 
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amount to 10. This may extend the disc spectrum to higher 
energies. 

Inside the disc spherisation radius, the disc wind sup- 
plies the funnel with gas starting from very deep regions 
JShakura fc Svunyaevl 19731; lEggum. Coroniti. fc Katj|l988l ; 



( snakura & bvunva evll^Aliiiggum. (Joromti. & Kataimsai; 
lOkuda et al.1 12003 ; lOhsuga et al.ll2005h . this gas is directly 



accelerated by the funnel radiation. The interaction of the 
gas accelerated to the velocities of ~ 0.2 — 0.3 c with new 
portions of the walls' wind creates the conditions quite suit- 
able for Comptonization. These conditions might be about 
the same throughout the whole funnel length where the jet 
is formed, as the high-speed outflow has about the same 
velocity. 

In the INTEGRAL observations a hard X-ray emis- 
sion (20 - 100 keV) with th e total hard luminosity o f 
10 35 erg/s has been detected ( Chcrcpashchuk et al. 2005), 



which may be well explained 1 Krivoshevev et al. 20091 ) 



terms of Comptonized thermal plasma surrounding the vis- 
ible jet base. 

Despite the fact that the jet formation and collima- 
tion mechanism is not clear yet, one may expect the Comp- 
tonization is significant inside the funnel. The observed spec- 
trum, taken when the di sc is the most open to the observer 
jKrivoshevev et all fcoOO). is nearly flat in the energy range 
of 10-20 keV. This spectrum may be a continuation of the 
flat 7-12 keV spectrum found here in the reflecting model. 

The soft excess detected in the SS433 spectra and in- 
troduced here as a BB component with the temperature of 
9 bb kO.1 keV (Tabled Fig.[H]l, may appear to be the direct 
radiation of the visible funnel wall at r > n (Fig. 3}. When 
we take into account the factor k « 3.2, which is a ratio of 
the total cone surface at r > r\ to the visible cone surface, we 
find the soft BB component luminosity L x ~ 3 ■ 10 37 erg/s. 

The soft BB component is shown in Fig. [8] Note that 
the luminosity of the soft component strongly depends on 
the adopted column density Nh ■ It is essential here that the 
soft c omponent is necessary dBrinkmann. Kotani. fc Kawail 
120051 . and herein) in order to understand the SS 433 spectra 
within the framework of the thermal jet model. 

In Fig.[S] we show two versions o f the funnel's non- 
absorbed spectrum in the MCF mod el (|Fabrika et al.ll2006l ; 
iFabrika. Abolmasov. fc Karpovll2007l ), at r > 1.74 • 10 11 
(the value of ri derived above) and r > 5 ■ 10 
ing these spectra, we adopted the wind's photosphere radius 
and the temperature of r p h = 10 12 cm, T p n = 50000 K re- 
spectively (|12p . and T(r) oc r _1 as it was discussed above. 
This illustrates that the MCF spectra might be responsible 
for the soft excess observed in SS 433. Nevertheless, one can 
say nothing else, exept that the additional soft component 
is present in SS433 spectrum. 

Th e soft excess is observed in t he spectra of 
ULXs Jstobbart. Roberts, fc Wilms! 120061 ; iRobertsI 120071 ; 
iPoutanen et al.l 20071 ) . their temperature is 6 ~ 0.1 keV 
and it contributes to the total ULX luminosity of ~ 10 %. 
The soft component luminosity in ULXs is about 10 38 — 
10 39 erg/s. If the ULXs are nearly face-on SS 433 stars, their 
intrinsic soft component luminosity has to be some less be- 
cause of the geometrical beaming of their supercritical discs. 
However, the soft X-ray luminosity is model dependent, it 
strongly depends on the adopted absorption value. 



cm 

10 cm. Deriv- 



7 CONCLUSION 

The goal of the present paper was to find indications of 
the supercritical accretion disc funnel in SS 433 X-ray ra- 
diation. The bolometric luminosity of the disc in SS 433 is 
~ 10 40 erg/s and all the luminosity must be initially released 
in X-rays. This is confirmed by the well-established kinetic 
luminosity of SS433 jets, which is ~ 10 39 erg/s, and the fact 
that the jets are formed in the deepest places of the disc 
funnel close to the black hole. The observed X-ray luminos- 
ity of SS433 is ~ 10 36 erg/s and it is the radiation of the 
cooling X-ray jets. Both the orientation of SS 433 and the 
visibility conditions do not allow us to see any deep areas of 
the funnel. However, the strong mass loss of the accretion 
disc gives us a hope for detecting some indications of the 
funnel radiation. 

We have analysed the XXM spectra of SS 433 with a 
well-known standard model of the adiabatically cooling X- 
ray jets, taking into account cooling by radiation. We con- 
firm that the thermal jet model reproduces the emission 
line fluxes quite well. When the disc is most open to the 
observer, the visible blue jet base is ro ~ 2 ■ 10 11 cm, the 
gas temperature at ro is 80 w 17 keV. The IS gas column 
density is Nh ~ 1.5 • 10 22 cm -2 , which is in good agree- 
ment with the IS extinction found in optical and U V ob- 
serva t ions. We confirm the previous find ing l|Kotani et al.l 
1 19961 ; iBrinkmann. Kotani. fc Kawail 120051) that the red jet 
is probably not seen (blocked) in the soft energy range of 
0.8 -2.0 keV. However, both the model with blocked the 
red jet portions, and the model with the whole red jet vis- 
ible (but with some higher Nh) give the same result, that 
the thermal jet model alone can not explain the soft contin- 
uum. We also confirm that Nickel is highly overabundant in 
the j ets (jKotani et al.l 1 19961 ; IBrinkmann. Kotani. fc Kawail 



tne i ei 
2005). 



We have adopted the luminosity Lk = 10 39 erg/s and 
the jet opening angle 2#j = 1.5° as well established and 
known. However, the derived jet visible base ro is too small, 
it has to be ~ 10 12 cm to satisfy the observed eclipses by 
the donor star. Using the scaling formula ro oc Mf/flj one 
may find that during these XMM observations the jet kinetic 
luminosity was Lk ~ 2 • 10 39 erg/s. This scaling should not 
change the thermal jet spectrum. 

We find that the thermal jet model alone can not repro- 
duce the continuum radiation in the XMM spectral range. 
Therefore we use the thermal jet model together with the 
REFLION ionized reflection model and find quite a good 
representation of the spectra. Introducing the reflection not 
only explains the continuum curvature at the energies of 
> 3 — 4 kev and the fluorescent line of quasi-neutral iron, 
it reproduces the broad absorption fea ture at ~ 8 keV re- 
cently detected bv lKubota et al] l|2007l) . We independently 
study three energy ranges (2-4, 4-7 and 7-12 keV) with 
the same thermal jet model to find the reflection parame- 
ters. We find that the ionization parameter (|10p is about 
the same in these ranges, £ ~ 300, which indicates a highly 
ionized reflection surface. The illuminating radiation photon 
index changes from flat, F m 2, in the 7-12 keV range to 
T ~ 1.6 in the range 4-7 keV and to T < 1 in the range 
2-4 keV. 

We conclude that the additional reflected spectrum is 
an indication of the funnel radiation. The illuminating radia- 
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tion spectrum is flat in the range of 7 - 12 keV. With multiple 
scatterings in the funnel the hard radiation may survive ab- 
sorption. The reflected luminosity in this 7-12 keV spectral 
range is L re /i ~ 10 36 erg/s. The softer (2 — 7 keV) part of the 
illuminating spectrum (Fig. [8} carries a trace of absorption. 
This is assumed to be due to the multiple scatterings in the 
funnel. It is important to note that an existence of He- and 
H-like absorption edges at about zero velocity is a manda- 
tory property of the funnel spectrum to be able to produce 
the observed jet velocity due to the line-locking mechanism. 

The supercritical accretion disc spectrum is expected 
to be flat, however, it extends up to a few keV only. Comp- 
tonization may extend and flatten the spectrum to higher 
energies. The interaction of the high velocity gas moving 
along the jet axis with the walls' wind makes the conditions 
for the Comptonization throughout the whole funnel length 
consistent. Direct evidences of the Comptonized radiation 
have recently been found i n the INTEGRAL observations 
l|Cherepashchuk et alj|2003 ). 

We have not found any evidences of the reflection in 
the soft 0.8-2.0 keV energy range. The soft excess is de- 
tected in the SS433 spectra. We suppose that in the soft 
X-rays we observe direct radiation of the visible funnel wall. 
We represented this component as a black body radiation 
with a temperature of On ~ 0.1 keV and luminosity of 
Lbb ~ 3 ■ 10 37 erg/s. The soft excess may be also fitted 
with a multicolour funnel (MCF) model. Both the soft X- 
ray luminosity and its derived spectrum strongly depend on 
the column density Nh- The soft component is necessary to 
add in order to understand the SS 433 spectra within the 
framework of the thermal jet model. 

The soft excess is observed in the spectra of ULXs with 
the soft component temperature of 8 ~ 0.1 keV, which is 
identical with that found in SS 433. If the ULXs or some 
of them are nearly face-on versions of SS 433, one may ad- 
just their soft X-ray componen ts to the outer funnel walls 
radiation (jPoutanen et al.lfeOQTl '). 
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